Exchange-based CNOT gates for singlet-triplet qubits with spin orbit interaction 
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We propose a scheme for implementing the CNOT gate over qubits encoded in a pair of electron 
spins in a double quantum dot. The scheme is based on exchange and spin orbit interactions and on 
local gradients in Zeeman fields. We find that the optimal device geometry for this implementation 
involves effective magnetic fields that are parallel to the symmetry axis of the spin orbit interaction. 
We show that the switching times for the CNOT gate can be as fast as a few nanoseconds for 
realistic parameter values in GaAs semiconductors. Guided by recent advances in surface codes, we 
also consider the perpendicular geometry. In this case, leakage errors due to spin orbit interaction 
occur but can be suppressed in strong magnetic fields. 
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I. INTRODUCTION 

Standard quantum computing [l[ is based on encod- 
ing, manipulating, and measuring quantum information 
encoded in the state of a collection of quantum two-level 
systems - qubits. Spin-1/2 is an ideal implementation 
of a qubit, since it is a natural two-level system, and 
every pure state of a spin-1/2 corresponds to a state of 
a qubit. For this reason, spins have been considered as 
carriers of quantum information in a variety of proposals 
0. The initial proposal Q called for spins in single- 
electron quantum dots electrically manipulated by the 
exchange interaction and local time-dependent Zeeman 
fields. A variety of other encoding schemes and manip- 
ulation techniques 04H1 rely upon encoded qubits. In 
these schemes, the simplicity of qubit states and minimal 
number of physical carriers of quantum information are 
traded for less stringent requirements for experimental 
implementations. On one hand, the alternative setups 
protect from the most common types of errors by de- 
coupling the computational degrees of freedom from the 
most common sources of noise, and therefore allow for 
longer gating times. On the other hand, in some alter- 
native setups the manipulation without fast switching of 
the local magnetic fields becomes possible. 

The optimization in the encoding and manipulation 
protocols is always guided by the state of the art in the 
experiments. Recent results suggest that spin qubits can 
reside in a variety of material hosts with novel properties. 
Quantum dots in graphene [l2j and carbon nanotubes 
[l3( are less susceptible to the decoherence due to nuclei 
and spin-orbit interaction. Spins in nanowires show very 
strong confinement in two spatial directions, and the gat- 
ing is comparably simple In hole systems, the car- 
riers have distinct symmetry properties, and show cou- 
pling to the nuclear spins of a novel kind [l4[ . Recently, 
the experiments in silicon [l5j have demonstrated coher- 
ent manipulation of spins similar to the one achieved in 
the GaAs-based nanostructures. Within these hosts, the 
manipulation techniques that use exchange interaction, 
spatially inhomogeneous time independent Zeeman split- 



ting, and nuclear hypcrfmc interactions are within reach. 
Despite these developments, GaAS remains a promising 
route to spin qubits due to the highly advanced experi- 
mental techniques developed for this material. 

Here, we study the implementation of the quantum 
ates on the encoded two-spin singlet-triplet (ST) qubits 
ITtl [ijj using resources that closely resem- 
ble the ones available in the current experimental se- 
tups. There, the application of the time-dependent elec- 
tric fields through the gates fabricated into a structure 
[IH are preferred to time-dependent local magnetic fields. 
In addition, the nuclear spins HEM! and inhomoge- 
neous magnetic fields [2(J Hlj are possible resources for 
spin control. In the setups based on semiconductors, the 
electrons or holes in the quantum dots are influenced by 
the spin orbit interaction (SOI), which can contribute to 
the control [22| . 

In this work we present a scheme for control of ST- 
qubits which uses switching of the exchange interaction 
as a primary resource. We consider the scheme that is 
optimized for the application of quantum gates in the net- 
work of quantum dots. The construction of the CNOT 
gate uses pulses of the exchange interaction as the only 
parameter that is time dependent. The exchange inter- 
action itself is not sufficient for the universal quantum 
computation over the ST-qubits due to its high symme- 
try. The additional symmetry breaking is provided by 
nonuniform, but static magnetic fields. These fields de- 
scribe the influence both of magnetic fields, provided by 
the nearby magnets, and of the coupling to the nuclear 
spins in the host material via hyperfine interaction. De- 
pending on the scheme used for the application of the 
quantum gates, the optimal geometry is either the one 
in which the magnetic fields point parallel to the axis of 
symmetry of the SOI or perpendicular to it. 

One major problem in the realization of two-qubit 
quantum gates (in particular, we consider here the 
CNOT gate based on conditional phase gates) , is the pos- 
sibility of leakage errors where the spin states defining the 
logical qubit leave the computational space. These errors 
move the state of four spins from the 4-dimcnsional com- 
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putational space of two qubits into some other portion of 
the 16-dimensional Hilbert space of four spin-1/2 parti- 
cles. We consider two ways of addressing this problem. 
One scheme possesses the axial symmetry due to the fact 
that the SOI vector and magnetic fields are parallel. For 
this 'parallel scheme' we are able to construct a perfect 
CNOT gate, if we are able to control all the available 
parameters. Having in mind 2D architectures, we also 
consider the CNOT gate between two qubits in the case 
when the SOI vector and magnetic fields are perpendic- 
ular to each other. Here we cannot prevent the leakage 
out of the computational space, however, we show that 
it is suppressed by a ratio between the SOI and Zeeman 
energy coming from a strong external magnetic field. 

All our constructions assume that the controlled inter- 
actions are switched in time by rectangular pulses. Any 
deviations from this form of time dependence lead to ad- 
ditional corrections and affect the fidelity of the gate. 

The paper is organzicd as follows. In Sec. [Hi we intro- 
duce the model for the double dots and effective Hamil- 
tonians for field gradients and exchange and spin orbit 
interactions. In Sec. IIII1 we consider the parallel geome- 
try and derive the CNOT gate via the conditional phase 
gates and swap gates, all based on exchange. There we 
also give estimates for the switching times. The scheme 
for the perpendicular geometry is then addressed in Sec. 
TVl and we conclude in Sec. [V] 



II. THE EFFECTIVE MODEL 

We consider in the following singlet-triplet (ST) qubits 
that are implemented by two electrons confined to a dou- 
ble quantum dot system [HQ, see Figs. [T]and [2] Such 
ST-qubits have been realized successfully in several labs 
(23j . and single and two-qubit operations have also been 




FIG. 1: Parallel geometry: four quantum dots (yellow discs) 
aligned along the x-axis in the presence of an external mag- 
netic field B that is applied parallel to the SOI vector j3 (red 
arrow), which must be perpendicular to the line of the dots 
and which we take to be the z-direction. At each dot, there 
is a local magnetic field Bi (blue arrows), also assumed to 
be parallel to B, but with alternating orientations as indi- 
cated. The direction of B defines the spin quantization axis. 
The dots are defined electrostatically by metallic gates (light 
green structures). Each dot contains a spin-1/2, and the ex- 
change (Jij) and the SOI-induced (/3y) interactions between 
the spins can be controlled by changing the electrostatic po- 
tential between the dots. The dots 1 and 2 (from left to right) 
define the first ST-qubit, and the dots 3 and 4 the second ST- 
qubit. 



demonstrated recently @, QU H3, 113 ■ There are sev- 
eral schemes for the fundamental CNOT gate, which can 
be divided into two classes, schemes which make use of 
exchange interaction and schemes which do not, but in- 
stead rely on coupling of dipole moments [lfj. The lat- 
ter schemes has the disadvantage to be rather slow and 
also to be affected by charge noise rather strongly. Here, 
we focus on exchange-based schemes specifically adapted 
to quantum dots in III-V semiconducting materials, that 
have the SOI, such as in GaAs or InAs quantum dots. Al- 
though the SOI is typically small compared to the level 
spacing of the dots, it needs to be taken into account in 
order to achieve high fidelities in gate operations. 

We focus now on two such ST-qubits and assume them 
to be encoded in four quantum dots that are arranged in 
a row, see Figs. Q]and [U The external magnetic field B is 
assumed to give the largest energy scale and determines 
the spin quantization axis z. The Hilbert space of four 
spins-1/2 is spanned by 2 4 = 16 basis states. The total 
spin of the system is given by S = y ■_ 1 Si, where is a 
spin-1/2 operator acting on the spin in a dot i = 1, 2, 3, 4. 
Due to axial symmetry, the z-component, S z , becomes a 
good quantum number, and the eigenstates correspond- 
ing to S z = span a six-dimensional subspacc. We define 
the computational basis of the two ST-qubits in this sub- 
space as 

100) Hint), in) =\nn), 

101) Him), |io) =imt), (i) 

where {0, l}® 2 denotes the ST-qubit space, 'f and 
denote states of the quantum dot spins corresponding to 
the projection Sf = ±1/2 on the z axis. The remaining 
two states, 

= ittu) . \h) = \im) , (2) 

belong to the non-computational leakage space. 

Besides the externally applied magnetic field B, we 
allow also for local magnetic fields B^ that are constant in 
time (at least over the switching time of the gate). Such 
local fields can be generated e.g. by nearby micromagnets 
[2ll | or by the hyperfine field [25M30I] produced by the 
nuclear spins of the host material. In the latter case, in 
order to reach a high fidelity, it is important to perform a 
nuclear state narrowing [25[ , i. e. to suppress the natural 
variance SBi ~ A/y/N ~ lOmT to a smaller value. In the 
ideal case, one should try to reach a limit where |B;| ~ 
50mT <C |B| and the fluctuations in B; are much smaller 
than |B,| . 

The magnetic fields B^ are assumed to point along the 
external field B, so that they preserve the axial symme- 
try of the problem. However, the Bj's should have differ- 
ent values (to create field gradients between the dots), a 
simple scenario being local fields of opposite directions on 
neighboring dots, see Figs. [I] and [2] The corresponding 
Zeeman Hamiltonian is given by 



(3) 
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FIG. 2: Perpendicular geometry: similar to the setup shown 
in Fig. [T]but with the difference that here the SOI vector /3 
(red arrow) is perpendicular to the magnetic fields B and B^, 
which define the spin quantization axis z. 



where the effective magnetic fields are defined in terms of 
energy as b = gfiBB and 6j = gptsBi, respectively, with 
g being the electron g-factor and [1b the Bohr magneton. 

The exchange interaction Jij(t) couples the electron 
spins of nearest neighbor dots i and j and can be con- 
trolled electrostatically @, HH HH . If the tunnel barrier 
between the dots is high we can treat them as indepen- 
dent. If the tunnel barrier is lowered and/or a detuning 
between the dots is applied, the two spins interact with 
each other, leading to an effective description in terms of 
an Heisenberg Hamiltonian 0, [22|, [3l[ , 

H ex = J 12 (t)Sx ■ S 2 + J 23 (t)S 2 ■ S 3 + Mt)S 3 ■ S 4 • (4) 

We assume that the magnetic field is sufficiently large 
compared to exchange energies, i.e. Jy <C B, to avoid 
admixture of triplets via the SOI (see the discussion of 
perpendicular geometry below). 

We note that ideally it is best to switch the exchange 
Jij by changing the corresponding intcr-dot barrier- 
height or distance, instead of detuning the double dot 
by a bias e. Detuning is harmful for two reasons. First, 
detuning can admix other unwanted states (for example, 
(0,2)S, see Fig. 2 in H2|). To analyze the errors coming 
from detuning, one needs to go beyond the effective spin 
Hamiltonian and consider the microscopic model for the 
double dots which includes Rashba and Drcsselhaus SOI 



and inhomogeneous fields, see Ref. |22|. Second, the con- 
trol of via the tunnel barrier preserves the symmetry 
of the charge distribution in the double dot and thus, 
in particular, avoids the creation of dipole moments. In 
contrast, such dipole moments are unavoidable for detun- 
ing, and in the regime with dJij/de ^ charge noise can 
enter most efficiently the qubit space, causing gate errors 
and decoherence of the ST-qubits [25|. Thus, symmetric 
exchange switching is expected to be more favourable for 
achieving high gate fidelities. 

Next, we account for the effects of spin orbit inter- 
action. Following Refs. [321434} we model the SOI by 
a Dzyaloshinskii-Moriya (DM) term for two neighboring 
quantum dots (see e.g. Eq. (1) in [HI), 



H™ 1 =f3 ij (t)-(S i xS j ), 



(5) 



where the SOI vector (t) is perpendicular to the line 
connecting the dots. First we consider a 'parallel geome- 
try' (sec Fig. [T]) where the SOI vectors (3^ are all parallel 



to each other and the magnetic fields B and B 3 are as- 
sumed to be parallel to the SOI vectors. This preserves 
the axial symmetry of the spin system, and by definition, 
we choose the direction of B to be the z-axis. This leads 
to 



H 



SOI 



(6) 



where the summation runs over neighboring dots i and 
j. The strength of the SOI, fy, depends on the distance 
between the dots as well as on the tunnel coupling be- 
tween them. This allows us to assume that both Jjj(t) 
and Pij (t) are switched on and off simultaneously [35l437| . 

We note here that both H ex and , being axially 

symmetric interactions, preserve the ^-component of the 
total spin S z . This means that our set-up is protected 
from leakage to the subspace with S z ^ 0. However, it is 
not protected from the leakage to the non-computational 
space given by Eq. @. By a proper design of gates this 
leakage can be minimized. 

Alternatively, in a 'perpendicular geometry' (see Fig. 
[5} the axis of the quantum dots is aligned parallel to 
the applied magnetic field, in the z-direction. The SOI 
vector [3^ is perpendicular to this, and we take it to be 
in x-direction. The corresponding Hamiltonian becomes 



Hf_ OI = J2Pij(t)(s?s?-sisy). 



(7) 



Here, the SOI vector (3^ breaks the axial symmetry of 
the system, and the total spin S z is no longer a good 
quantum number. As a consequence, leakage into the 
non-computational space S z ^ is possible, in principle. 
However, this coupling involves higher-energy states and 
can thus be suppressed by choosing a sufficiently large 
magnetic field such that j3/B <C 1. In contrast, the SOI 
does not couple states within the computational space, 
since the matrix elements of Hf 01 between the states 
||t) and WX) vanish. 



III. PARALLEL GEOMETRY 

In this section we concentrate on the parallel geometry, 
see Fig. [TJ Using the axial symmetry of the problem we 
are able to construct a sequence of gate operations that 
implements the CNOT gate [l|, defined on the logical 
ST-qubits given in Eq. ([T]) by UcNor\a,b) = \a, a ©6), 
where a, b = 0, 1. 

One important step in this construction is the imple- 
mentation of the 7r/4-gate U^/i (see discussion below and 
[H). For this gate we propose the following scheme con- 
sisting of four steps, 



C*23 — > (^12, ^34) -> C23 — > (7ri2, ^34), 



(8) 



where the conditional phase gate C23 exchange-couples 
the dots 2 and 3 and adds a phase factor to the two ST- 
qubits (see below). The swap gates 7Ti2 and ^34 exchange 
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spin states on the dots 1 and 2 and the dots 3 and 4, and 
can be performed in parallel. 

A major issue in the implementation of U^/a is the 
avoidance of leakage errors during the coupling of qubits. 
To keep qubits in the computational space, operations on 
spins 2 and 3 must be constructed in such a way that the 
resulting gate is diagonal in the basis given by Eq. (JTJ. 
This can be achieved in two ways. The first approach is 
to use strong pulses that lead to fast rotations around the 
Bloch sphere. The second approach is to use adiabatic 
pulses that are protected from the leakage to states with 
different energies f38| . However, the adiabaticity requires 
a longer pulse time. In the present work we focus on the 
first approach. 



A. Conditional phase gate C23 




FIG. 3: The Bloch sphere is defined with the north pole cor- 
responding to \\\) and the south pole corresponding to \\."\). 
The effect of the unitary evolution operator Ug [see Eq. (|17l) ] 
on a state in the space {|t^)i 14-T)} is equivalent to the rotation 
on the Bloch sphere around the vector J23e x — P2ze y + Af>23e- . 
The conditional phase gate C23 corresponds to the full rota- 
tion on the Bloch sphere (shown by blue and green circles). 



In this subsection we describe the phase gate, C23, in- 
volving the exchange and SOI interactions only between 
dot 2 and 3, while dots 1 and 4 are decoupled from dots 
2 and 3, i.e., J\ 2 = J34 = and /?i2 = (3 3 a = 0. In this 
case, the effective Hamiltonian is given by 



H c = H B + H ex 



Hi 01 = H° 



H 



23 



(9) 



where we present it in block-diagonal form. The part 
of the Hamiltonian Hg = (b + bi)S* acts only on spins 
located at the dot i = 1,4. The other part of the Hamil- 
tonian, on spins located at the dots 2 and 3, 



Hg = (b + b 2 )S z 2 + (b + b 3 )SI + J23S2 • S 3 



(10) 



Here, we assume a rectangular pulse shape for the ex- 
change and spin orbit interactions, and from now on we 
treat J23 and j3 23 as time-independent parameters. In 
this case, the unitary gate U c is a simple exponential of 
the Hamiltonian, 

jjC = e -iH c T c = e -iH°Tc e -iH&Ta e -iH°T a> ( n ) 

The spins of dot 1 and 4 do not change in time apart 
from a phase factor coming from the corresponding mag- 
netic field. In contrast, the spins of dot 2 and 3 change in 
time and acquire phases, as we describe next. For this we 
express Hg 3 as matrix in the basis {|tt)j lit) > 




Hg 

hSl 



(12) 



where the block-diagonal form reflects the conservation 
of S z . In the case of two parallel spins, the correspond- 
ing Hamiltonian is given by Hg = J23/4 ± (62 + b 3 )/2, 
which just assigns a phase to the spins. In the case of 
antiparallel spins, the Hamiltonian Hg is given by 



Hg = ^ 



1 /-J 23 /2 + A6 ; 



23 



2.3 



^23 + ^23 

-J23/2- Ab 23 



(13) 



where Ab 23 = b 2 — b 3 . Here, Hg describes the cou- 
pling between the states |fi) and |4t)i which, in general, 
leads to leakage errors. This leakage can be prevented by 
choosing the pulse duration, Tq, in such a way that the 
corresponding unitary operator Ug = exp[— iHg Tc] is 
diagonal in the basis {|ti)j lit)}- It is more convenient to 
consider the evolution given by Hg on the Bloch sphere. 
For that we rewrite Hg in terms of pseudospins, 



Hg = -J23/4+(J 23 /2) n 23 -T, 

^23 = V^+Z^ + ^M 2 , 
»23 = (^23, — /?23, Ab 2 gj / J 23 , 



(14) 
(15) 
(16) 



where the unit vector 1123 defines the rotation axis on the 
Bloch sphere, see Fig. [31 and the pseudospin, acting on 
the states {\ti), lit)}, is described by the Pauli matrices 
r. The north pole corresponds to |tl) and the south pole 
to The exchange interaction, J23, being the largest 
scale in Hg , forces n 2 3 to be aligned mostly along the 
x-axis. If we neglect the SOI and any field gradients, the 
rotation on the Bloch sphere takes place in the j/z-planc. 
In the presence of SOI and field gradients the rotation 
axis 1123 deviates from the x-axis. 

The unitary time evolution operator Ug corresponding 
to Hg takes the form 

Ug = exp(iac / '2) (cos ac + i n 23 ■ Ts'ma c ), (17) 

where ac = J2 3 Tc/2. The duration of a pulse is deter- 
mined by the condition that we obtain full rotations on 
the Bloch sphere (see Fig. [3]), 



Tr 



2ttN c 
J: 



(18) 



23 



where Nc is a positive integer. Note that deviations from 
Eq. (fT8"|) lead, again, to leakage errors. 

As a result, the qubit states with parallel spins on dot 
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2 and 3, acquire the phases 

0oi = ^(J 23 /2 -b 1 + b 2 + b 3 - b 4 )T c , 

010 - ^(J 23 /2 + h-b 2 - b 3 + b 4 )T c , (19) 

while the qubit states with antiparallcl spins on dot 2 
and 3 acquire the phases 

011 = \{bx - b 4 - (-l) Nc J 23 /2)T C , 

0oo = \{~bi + b 4 - {-\) N °J 23 /2)T C . (20) 

Here, a b corresponds to a phase acquired by a two-qubit 
state \ab) . We note that phases produced by the magnetic 
fields terms will be canceled during the second C 23 pulse 
after the 7r-pulses have been applied to the qubits. 

B. Swap gates ni2 and 7134 

In this subsection we discuss the swap gates n\ 2 and 
tt 34 that exchange spin states between dot 1 and 2, and 
dot 3 and 4, respectively. The swap operation is a one- 
qubit-operation, so dots 2 and 3 should be decoupled 
during the swap pulse, i.e., J 23 = and /3 23 = 0. The 
swap gate ir = (7112,^34) is implemented, again, by a 
rectangular pulse and all parameters are assumed to stay 
constant during the switching process. This simplifies 
the unitary evolution operator, 

W = U? 2 XJl A = e - lH ^e- lH ^. (21) 

Further, we focus on the first ST-qubit (dots 1 and 2) and 
consider only tt 42 (ir 34 can be obtained analogously). We 
also note that since S z is conserved, the final state is al- 
ways given by a linear combination of the states |4_t) and 
|t4-) ° n dot 1 and 2. Within this subspace the effective 
Hamiltonian is given by 

_ 1 f-Ju/2 + A6 12 Ji 2 +if3 12 \ ( 
12 ~ 2 V J12 ~ ijSia -J12/2 - Ab 12 ) ' W 

or in pseudospin representation [see Eq. (|14p] , 

Hl 2 = -J 12 /4 + (J 12 /2) n 12 ■ r, (23) 

Jia = /M + W, (24) 

nu = (Ji2,-/3i2,A6i 2 )/Ji2, (25) 

where A&i 2 = b\ — b 2 and the Pauli matrix r, acts in 
the pseudospin space spanned by \\\) and Again, 
the unit vector n 12 defines the rotation axis. The corre- 
sponding unitary evolution operator reduces to the form 

U12 = exp(ia w /2)(cosa T + mi 2 • rsina T ), (26) 

where a, = 




FIG. 4: Two schemes for the swap gate, tti 2 : |t4-) ~~ Hit): 
composed of three consecutive rotations on the Bloch sphere 
where the offset induced by A612 and P12 is fully compen- 
sated, a) First, starting from the north pole, we turn on J12 
and /3i2 to induce rotation around Ji2e x —/3i2e y + A6i 2 e z (up- 
per blue arc) until we reach the equator where we turn off J12 
and /3i2- Second, we let the state precess around the z-axis 
in the equatorial plane until the mirror point of the starting 
point on the equator is reached (brown arc). Third, we induce 
once more rotation around Ji2e x — Pi2e y + Ab^e^ (lower blue 
arc) until we reach the south pole. Lower panel: associated 
rectangular gate pulses and switching times j^ al ' a2 ' a3 ' f or the 
three steps, b) Alternative scheme where during the second 
step the state precesses along the equator until it reaches the 
point diametrically opposite to its starting point (brown arc). 



The swap operation should exchange the states |tl) 
and \it)- 111 t ne absence of SOI and field gradients, the 
unitary evolution operator U* 2 corresponds to a rotation 
around the a;-axis (ni 2 = e x ) in the yz-plane. At half 
the period, T° = 7r/2Ji 2 , a state evolves from the north 
to the south pole and vice versa, z.e., U 42 \t® oc t x . 

However, in the presence of SOI and/or field gradients, 
the rotation axis ni 2 is not aligned with e x (compare with 
Fig. [3]), so the trajectory starting at the north (south) 
pole would never go exactly through the south (north) 
pole. The corresponding deviations lead to errors in the 
7Ti 2 gate on the order of /Sb\ 2 + /3 2 2 /Ji 2 . This means 
that it is impossible to produce a perfect swap operation 
with only one single rectangular pulse. However, by ap- 
plying a sequence of several rectangular pulses, it is still 
possible to produce a perfect 7Ti 2 gate, as we demonstrate 
next. 

Indeed, this goal can be achieved by the following three 
steps (see also Fig. [4}. First, we switch on the exchange 
interaction J 12 between the dots 1 and 2 (this also au- 
tomatically switches on the SOI /3i 2 ). We assume J\ 2 to 
be larger than A&i 2 and/or /3i 2 . The rotation axis \\\ 2 in 
polar coordinates is given by 

ni2 = (sin ■& cos <p, sin 1? sin ip, cos (27) 

From now on we work in the coordinate system in which 
the x-axis points along J\ 2 e x — /?i 2 e J/ , so tp — and 
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< cos-d = A612/J12 < l/y/2. The strength and the 
duration of the rectangular pulse is chosen in such a way 
that the rotation reaches the equator of the Bloch sphere. 
The initial and final vectors on the Bloch sphere are given 
by 

& = (0,0,1), (28) 
(t> } = (cottf, v / sin(2tf-7r/2)/siritf,0), (29) 

allowing us to find the rotation angle al a1 ^ = 7r/2 + 
arcsin(cot 2 and the corresponding pulse duration 

T^ Ql) = [tt + 2 arcsin(cot 2 •ff)]/2J 12 . (30) 

Second, after switching off exchange and spin orbit in- 
teractions, J12 = and (3\2 = 0, the rotation takes place 
around the z-axis in the equatorial plane, at a precession 
frequency determined by the field gradient A&12. The 
rotation angle becomes or^ = 2[tt — arccos(cot 1?)], and 
the pulse duration is given by 

T W =a W/Ab 12 . (31) 

Finally, we repeat the first step by applying a pulse of 
the same strength J12 (P12) and during the same time, 

rp(a3) rp(al) 

An alternative scheme (b) is presented in Fig. 2}d. 
During the second step the state evolves on the Bloch 
sphere only over half of the equator, a^ 2 ^ = n, with the 
corresponding pulse duration T^ 2 ^ = oi^ 2 " 1 j ISbyi- The 



duration of the third pulse is given by 

Tf 3) = [3tt - 2 arcsin(cot 2 0)]/2J 12 . (32) 

The second step is the slowest one in these schemes, so 
the scheme (b) has an advantage over the scheme (a) by 
being faster as it requires less rotation on the equator. 
However, scheme (b) requires better control of parame- 
ters, since Ti bl) ^Ti b3) . 

Here we note that it is also possible to switch off the 
exchange coupling J12 not only on the equator but also 
at any other point on the Bloch sphere. Moving away 
from the equator speeds up the gate performance, but 
demands greater precision in the tuning, since the rota- 
tion proceeds along a smaller arc and in shorter time. 

The scheme presented above confirms that it is possible 
to construct a perfect 7r swap gate even in the presence of 
z-component of the SOI vector, /3 12 , or local field gradi- 
ents, Abi2, by adjusting the pulse durations. The other 
two components of the SOI vector couple states of differ- 
ent total spin S z , and thus cause leakage errors. There- 
fore, it is optimal to orient the magnetic fields (defining 
the spin quantization axis z) along the SOI vector, /3 12 . 



C. CNOT gate 

After the execution of the four-step sequence given by 
Eq. (JSJ an initial qubit state is restored but with a phase 
factor [see Eqs. (fTS) and ([2"D])]: 



[00] e^°°[(K)] ( " 12 ' " 34) > e ^°°[ll] ^ e ^+M [n] ^\ [ 0]e^ 11+ *»\ 

[11] e *"[H] (7ri2 ' e **"[00] ^> e ^ 11+< «[00] ( " 12 ' " 34) > [ll] e *<*"+*»>, 

[01] ^ e ^[01] ( " 12 ' " 34) > e^[ll] ^ e *(^+^o)[H] ( " 12 ' [QQ]e i ^+^\ 

[10] e*°[10] ( " 12 ' " 34) > e*°[01] ^ e <(*>i+*w)[ 01 ] " 34) > [ 10 ] e *(«oi+*io) ) (33) 



with 0oi + 0io = -(0n + 0oo) = J23TC/2, where we 
assume Nc [see Eq. (|18l) ] to be even. The total gate 
acting on the qubits as defined by Eq. (|33j) can be written 
in the compact form 

e -i(J 23 T c /2)crla^ 

where the Pauli matrix a\ acts on the first ST-qubit 
(formed by dot 1 and 2) and a\ on the second one (formed 
by dot 3 and 4) with j = x,y,z. Choosing 

Ib= (47rm- |)/J 23 , (35) 

where m is a positive integer, we obtain the 7r/4-gate, 

U v/A = e i tt'*. (36) 



Both Eqs. (fTS)) and (|35l) should be satisfied simultane- 
ously For example, if m = 1 and Nc = 2, we get 

^23 = ^V^+W 3 - (38) 

From this we can estimate the total switching time to per- 
form the 7r/4-gate that is a sum of the switching times 
at each step [see Eq. flSJ)]. For the scheme discussed in 
Subsecs. IIII Al and IIII Bl the slowest part is given by the 
swap gates 7ri2 and W34, whose switching time is limited 
by field gradients (due to nuclear spins SEMI and/or 
micromagnets 0,[ll|), |A6 12 | = |A6 34 | « 10 mT, which 
corresponds to T^ 62 ' w 10 ns. The gate can be faster 
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./ 23 (t) 



FIG. 5: An alternative scheme for the conditional phase gate 
C23 (see also Fig. [3]). Instead of the rotation with the pulse 
defined by Eqs. (|18|) and (|35p we consider a sequence of three 
pulses. During the first and the third pulses (blue arcs) the 
state precesses quickly acquiring the 7r/4-phase, see Eq. (|35|) . 
During the second pulse, J23 = 0, the state precesses around 
the z-axis over a shorter path than the one in Fig. [3] As a 
result, the switching is faster. 



if the rotations around the z-axis are performed not on 
the equator but more closely to the poles. This allows 
us to decrease the switching time of the swap gate to 
2 ns; however, this would require a more precise con- 
trol over the pulses. The same trick can be used to de- 
crease the switching time of the conditional phase gate 
C23 (compare Figs. [3] and [5]). If field gradients larger 
than 10 mT are used, the switching rates will be propor- 
tionately larger. 

Using the 7r/4-gate, we construct the controlled phase 
flip (CPF) gate Ucpf = diag(l, 1, 1, — 1) (see footnote 
[13] in 0) as 



U C PF = U 7v/A e- i ^+^- 1 l 
Finally, we obtain the CNOT gate, 

I ' 



(39) 



u C NOT = r , (40) 



IV. PERPENDICULAR GEOMETRY 

In the previous section we have discussed the parallel 
geometry for which we were able to construct a perfect 
CNOT gate under the assumption that we have a com- 
plete control over the parameters. The CNOT gate, to- 
gether with singlc-qubit gates, allows us to simulate any 
other quantum gate and its implementation is a crucial 
step toward the realization of a quantum computer [l[. 
In a next step, many such elementary gates need to be 
connected into a large network. In recent years, the sur- 
face code (39j has emerged as one of the most promising 
platforms for this goal due to its large threshold of about 
1% for fault tolerant error correction [HI. This platform 
requires a two-dimensional (2D) geometry and can be im- 
plemented in semiconductors of the type considered here. 

The basic 2D scheme is illustrated in Fig. [6l There, 
we show an array of quantum dots where two neighbor- 
ing dots in a given row represent one ST-qubit. These 
quantum dots are embedded in a semiconductor where 
the SOI is of the same type in the entire structure. As a 
typical example we mention Rashba and/or Dresselhaus 
SOI that both depend on the momentum of the electron. 
As a result, the direction of the SOI vector f3 is always 
perpendicular to the line along which two quantum dots 
are coupled by exchange interaction [22|. Thus, for the 
coupling of two such qubits in the same row, the SOI 
vector (3 is parallel to the external magnetic field B, so 
we can use the scheme designed for the parallel geometry 
in the previous section. In contrast, if we want to couple 
two qubits from neighboring rows, we should also con- 
sider a 'perpendicular geometry' (see Figs. [T] and H>]) in 
which the SOI vector (3 is perpendicular to the magnetic 
field B. 

This perpendicular geometry is characterized by sev- 
eral features. As was mentioned before, the axial sym- 
metry in this case is broken by the SOI, leading to 
the coupling between computational (S z = 0) and non- 



by using the CPF gate and performing a basis rotation on 
qubit 2 (a single-qubit rotation by 7r/2 about the y-axis), 



Ucnot = e l * a i U C PFe~ 



(41) 



In summary, the full sequence of operations for the 
CNOT gate, Ucnot, is given by 



e 1 *" 2 [(7ri27r 3 4)C23(7ri27r34)C'23]e 



-if (ctJ+o-J) p-tf a\ 



(42) 



We note again that this result has been derived un- 
der the assumption of rectangular pulse shapes. This is 
certainly an idealization, and in practice we expect devi- 
ations from this shape to cause errors for the gates and 
to affect the gate fidelity The study of this issue, being 
very important for practical purposes, requires a separate 
investigation and is beyond the scope of this work. 



B 




I I I . ; I 3 . j 



FIG. 6: Schematic setup for 2D architecture. Two dots define 
an ST-qubit (black dotted ellipses). An external magnetic 
field B and local magnetic fields (blue arrows) are parallel. 
In the case of the coupling between two qubits from the same 
row (red solid ellipse) the SOI vector /3|| is parallel to the 
magnetic field B, corresponding to the parallel geometry (see 
Fig. [1} . In the case of coupling between two qubits from two 
neighboring rows (red dashed ellipse) the SOI vector f3 ± is 
perpendicular to the magnetic field B, corresponding to the 
perpendicular geometry (see Fig. [2]). 
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computational (S z 7^ 0) space. If the magnetic field B 
is sufficiently large to split the triplet levels T± far away 
from the computational space (g/iBB ^s> /?), we can ne- 
glect this leakage. We estimate for GaAs B = 5T rj 
100/^eV/g/iB- At the same time, the SOI does not couple 
states within the computational space, so for the realiza- 
tion of the phase gate C23 and the swaps gates 7Ti2 and 
7T34 we can use the same scheme as in Sec. Illll onlv with 
= 0. 

V. CONCLUSIONS 

We have studied the implementation of the CNOT gate 
for ST-qubits in a model that is appropriate for current 
experiments @, [Ili H3| ■ The setup consists of an array 
of quantum dots with controlled growth direction and 
the relative orientation of the dots. Pairs of neighboring 
dots build the ST-qubits, where the quantization axis 
is determined by an externally applied magnetic field 
B. Moreover, we introduce an inhomogeneity in mag- 
netic fields, Bi, by local micromagnets or by the hyper- 
fine field. The resources used for time-dependent control 
are the exchange interaction Jij(t), the field gradients 
AB i:j (t), and the SOI vector (5 t] {t). 

If the SOI vector (3 is parallel to the external (B) and 
local magnetic fields (B 2 ), we are able to construct a per- 
fect scheme for the CNOT gate based on the 7r/4-phasc 
gate, U„u, consisting of four basic steps. Two of the 
steps involve interaction of spins that belong to different 



qubits, and open the possibility of leakage errors. Un- 
der condition of total control over system parameters, 
we show that the leakage can be eliminated. In the other 
two steps, the tuning of exchange interaction enables per- 
fect swap gates even in the presence of field gradients and 
SOI. 

Motivated by recent results on the surface code, we 
shortly comment also on the 2D architecture. Here we 
encounter a situation in which the SOI vector (3 and the 
magnetic fields, B and B^, arc perpendicular. In this 
case, the leakage to the non-computational space with 
S z 7^ is inevitable. However, it can be made very small 
as long as fi/gfisB <IC 1. 

Depending on the pulsing scheme, the switching times 
for the conditional phase gate are shown to lie in the 
range 1-100 ns for typical GaAs parameters. Compared 
to the experimentally established decoherence times of 
about 250 /j,s [l7|, this switching is sufficiently fast and 
shows that a CNOT gate based on exchange is a promis- 
ing candidate for experimental realizations. 
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